BNL-NT-06/14 
RBRC-574 

Single Transverse-Spin Asymmetry in Drell-Yan Production 
at Large and Moderate Transverse Momentum 

Xiangdong Ji/'|j Jian-Wei Qiu,^'[] Werner Vogelsang,^' ^'|f| and Feng Yuan^'lj 

^Physics Department, University of Maryland, College Park, MD 20742 
^Department of Physics and Astronomy, 
. . . , Iowa State University, Ames, I A 50011 

O ■ ^Physics Department, Brookhaven National Laboratory, Upton, NY 11973 

CN ; ^RIKEN BNL Research Center, Bmldmg 510 A, 

Qh' Brookhaven National Laboratory, Upton, NY 11973 



< 



O 

o 

o 



(Dated: February 2, 2008) 

Abstract 
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f^ ■ We study the single-transverse spin asymmetry for the Drell-Yan process. We consider pro- 



duction of the lepton pair at large transverse momentum, q± ~ Q, where Q is the pair's mass. 
The spin asymmetry is then of higher twist and may be generated by twist-three quark-gluon 
correlation functions. Expanding the result for q± <C Q, we make contact with the transverse- 



^j^. momentum-dependent QCD factorization involving the so-called Sivers functions. We find that 

p. 



QJ , the two mechanisms, quark-gluon correlations on one hand and the Sivers effect on the other, 



contain the same physics in this region. This ties the two together and imposes an important 
^ ■ constraint on phenomenological studies of single transverse spin asymmetries. 
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I. INTRODUCTION 

Single-transverse spin asymmetries (SSAs) play an important role for our understanding 

of QCD and of nucleon structure. They have a long history, starting from the 1970s and 

1980s when surprisingly large SSAs were observed in hadronic reactions such as p^p — > vrX 

n 

at forward angles of the produced pion 1]. The last few years have seen a renaissance in the 
experimental studies of SSAs. The HERMES collaboration at DESY, SMC and COMPASS 
at CERN, and the CLAS collaboration at the Jefferson Laboratory have investigated SSAs 
in semi- inclusive hadron production e A^^ -^ enX in deep- inelastic scattering [^ . For proton 
targets, remarkably large asymmetries were found. The advent of the first polarized proton- 
proton collider, RHIC, has opened new possibilities for extending the studies of SSAs in 
hadronic scattering into a regime where the use of QCD perturbation theory in the analysis 
of the data appears to be justified. The STAR, PHENIX and BRAHMS collaborations 
have presented data for single-inclusive hadron production, and large single-spin effects at 
forward rapidities were found to persist to RHIC energies [2|. 

The observed large size of SSAs in hadronic scattering has presented a challenge for 



QCD theorists 



'wo mechanisms have been proposed [5|, la, ll| and been extensively 
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applied [7|, |8|, |9|, |lO|, luj in phenomenological studies. The first relies on the use of transverse- 
momentum dependent parton distributions for the transversely polarized proton. For these 
distributions, known as Sivers functions |^] , the parton transverse momentum is assumed to 
be correlated with the proton spin vector, so that spin asymmetries naturally arise from the 
directional preference expressed by that correlation. The other mechanism (referred to as 
Efremov-Teryaev-Qiu-Sterman (ETQS) mechanism) is formulated in terms of the collinear 
factorization approach and twist-three transverse-spin-dependent quark-gluon correlation 
functions of the proton p, |7|. 

A concept common to both mechanisms is the factorization of the spin-dependent cross 
section into functions describing the distributions of quarks and gluons in the polarized 
proton, and partonic hard-scattering cross sections, calculated in QCD perturbation theory. 
The question of which mechanism should be used in the analysis of a single-spin asymmetry 
is therefore primarily tied to the factorization theorem that applies for the single-spin ob- 
servable under consideration. For example, for the single-inclusive process p^p — > vrX, there 
is only one hard scale, the transverse momentum pt of the produced pion, and the SSA is 



power-suppressed ("higher-twist") by 1/pt- In this case, one can prove a co//mear factoriza- 
tion theorem in terms of the quark-gluon correlation functions [7[ , and the ETQS mechanism 
apphes. On the other hand, the observables typically investigated in deep-inelastic lepton 
scattering (DIS) are characterized by a large scale Q (the virtuality of the DIS photon) and 
by the much smaller, and also measured, transverse momentum qx of the produced hadron. 
In this two-scale case, single-spin asymmetries may arise at leading twist, i.e., not suppressed 



by 1/Q. The relevant factorization t 
dependent {ryiV)) functions [l2, 
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leorem is formulated in terms of transverse-momentum- 

n 



IJ, |15| , in particular the Sivers functions. We note that 
much progress has been made recently in understanding the underlying theoretical issues 
in the TMD QCD factorization. For example, the gauge- invariance properties of the nec- 



17,hJ,UJ for DIS and Drell-Yan 



essary TMD parton distributions have been clarified [16|, 
)rocesses, and an extension to more complicated hadronic processes has been considered in 

From what we just discussed, the two mechanisms for single-spin asymmetries might 
appear to be essentially unrelated. However, one can make an argument that a consistent 
theoretical description of the SSA for a hard process over its full kinematical regime requires 
both mechanisms to be present and to contain the same physics in the region where they 
both apply. To give a specific example - which will in fact be the example treated in this 
paper - let us consider the SSA for the Drell-Yan process when the invariant mass Q of the 
pair as well as its transverse momentum g^ are measured. We stress that we consider cross 
sections differential just in Q"^ and g_L, but not specifically in the angular distributions of 
any of the individual leptons. 

At relatively large pair transverse momentum, g_|_ ~ Q, there is only one large scale, and 
the SSA will be power-suppressed in that scale. This directs us to use the ETQS mechanism 
with its collinear factorization involving the twist-three quark-gluon correlation functions 
and corresponding hard-scattering functions calculated at lowest order from partonic 3^2 
processes. As an aside, when taking the real-photon (Q^ -^ 0) limit, this calculation will 
yield the result for the SSA in direct-photon production [3|. 

We can next investigate what happens when we make the ratio q^jQ small, keeping 
both scales perturbative, Q ^ g_L ^ Aqcd- We refer to gj_ in this regime as "moderate" 
transverse momentum. The ETQS mechanism will still apply here (even though the hard- 
scattering functions will develop large logarithms of the ratio q±/Q that will eventually 



need to be resummed to all orders in the strong coupling). At t 



le same time, however, 
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IJ], which involves the 



the factorization in terms of TMD distributions applies now [12, 

Sivers functions. If both mechanisms are internally consistent, they must describe the same 

physics in this region. 

In a recent publication ^21], we have demonstrated that the two mechanisms indeed 
provide the same description of the single-spin asymmetry for the Drell-Yan process in the 
regime Aqcd ^ q'± ^ Q, and that there is a direct correspondence between the Sivers 
functions and the twist-three quark-gluon correlation functions, derived also earlier in |l9|. 
The key observation is that, at moderate transverse momentum, the Sivers function may be 
calculated perturbatively, using the twist-three quark-gluon correlation functions. In other 
words, the ETQS mechanism generates a no n- vanishing Sivers function in this kinematic 

m 




regime. Although there have been earlier efforts to link the two mechanisms 
clear connection, at the level of physical observables, has been established only in 
the present paper, we present details of the derivation of our results in [21 1. 

Our results may in some sense be viewed as establishing a unification of the two mech- 
anisms. At large q±, the ETQS mechanism applies. At moderate transverse momentum, a 
smooth transition from the ETQS mechanism to the one based on TMD factorization occurs, 
with the two approaches containing the same physics. At yet lower q± (~ Aqcd), the TMD 
factorization still applies, containing in a natural way the transition from perturbative to 
non-perturbative physics. We believe that this unified picture should prove to be the best 
approach to phenomenological studies of SSAs. 

The remainder of the paper is organized as follows. In Sec. II, we introduce the twist- 
three quark-gluon correlation function and explain its physical significance. In Sec. Ill, we 
calculate the single-spin differential cross section for Drell-Yan production using the ETQS 
mechanism, and we study its limit q± -^ Q in Sec. IV. A comparison with the TMD 
factorization approach is made in Sec. V, and we conclude the paper in Sec. VI. 

II. TRANSVERSE-SPIN-DEPENDENT QUARK-GLUON CORRELATION 

Consider a transversely polarized proton traveling at nearly the speed of light. Its internal 
color electric and magnetic fields then have preferred orientations in the plane transverse 
to the proton's direction of motion. By parity invariance, the color electric field must be 
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FIG. 1: A vertex diagram representing twist-three quark-gluon correlation function in a proton. 
When the proton is polarized, the gluon field has a definite polarization. 

orthogonal to the spin of the proton. If averaged over the proton wave function, the field 
vanishes because the proton is color- neutral (and also because of time- reversal symmetry). 
However, if one multiplies the color field with the quark color current, the average may 
be non-zero. This average defines a quark-gluon correlation function that characterizes 
a property of a polarized proton. In some sense, the correlation describes how strongly 
polarized the color electric field is in a spinning proton. 

In this paper, we are mostly interested in the so-called light-cone correlations. For these, 
quark and/or gluon fields are separated along the light-cone direction ^~ (if S,^ denotes 
a space-time coordinate, the light-cone variables are defined as ^^ = (^° ± ^'^)/-\/2, ^_l = 
(^^)'C^))- If is the light-cone correlations that characterize the structure of the proton in 
high-energy scattering. One finds the following expression for the spin-dependent quark- 
gluon correlation described above J7[: 
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visualized by the diagram in Fig. 1. Here, \PS) is the proton state with four-momentum 
P^ and transverse-polarization vector S^ = (0,0, Sj,), ip, ip are the quark fields with color 
indices i,j and Dirac indices p, a, and F^°' is the gluon field strength tensor with octet color 
index a. Xi = k^-^/P^, Xg = k^ /P^, and X2 = k'^2/ P^ = Xi + Xg are the fractions of the 
proton's light-cone momentum carried by the quarks and the gluon, where their momenta 
as shown in Fig. 1. Furthermore, C is the light-cone gauge link, 

K2 



m2,^i) = pew 



-ig I dCA+iC 



(2) 



which makes the correlation operator gauge-invariant. Note that our sign convention for the 
strong couphng constant g follows from D^ = 5^ + igA^ for the covariant derivative. Also 
note that we have included a factor g in the matrix element in Eq. (^, as compared to the 
definition in \1\. In the second line of (Q), we have expanded the matrix in Dirac space, 
keeping only the leading, twist-three, term that contributes to the SSA in the Drell-Yan 
process, and neglecting contributions of yet higher twist. Here, Nc = 3 is the number of 
colors, Cf = (A^c* ~ ^)/'^^c = 4/3, and T^" denote the generators of the SU(3) color gauge 
group, e" is the 2-dimensional Levi-Civita tensor with e^^ = 1. 

Inverting Eq. (1), we can express the correlation function Tf{xi,X2) as follows: 

TFix„X2) = |^^^^e-^^"^-e^(---)^"^-ef 5^^(P^|^(0)7+^F+.(r)^(r)|P5> , 

(3) 
where the sums over color and spin indices are now implicit and, for simplicity, we have 
omitted the gauge link which vanishes if one chooses to work in the light-cone gauge, A'^ = 0. 
Note that we have suppressed a renormalization scale dependence of Tf- Equation Q 
provides the normalization for our calculations in the following sections. Because of parity 
and time-reversal invariance, the correlation function has the symmetry property [3] 

Tf(xi,X2) = Tf{x2,Xi) , (4) 

which will be used later on to simplify our results. 

III. SINGLE-SPIN DRELL-YAN CROSS SECTION 

We shall now calculate the single-transverse-spin dependent cross section for the Drell- 
Yan process, using the ETQS mechanism. We will first perform the calculation of the cross 
section at parton level. Convolution with the twist-three correlation function introduced 
in the last section and with the usual Feynman parton distributions for the unpolarized 
proton will then give the physical hadronic cross section. We shall limit ourselves to a 
situation where the dominant contributions come from the twist-three correlation function 
for a valence quark, and from an anti-quark or gluon of the unpolarized proton. In this way, 
we do not need to worry, for example, about other types of twist-three correlation functions. 



such as purely gluonic ones 



2J|. In an actual experiment, such a situation may be realized 






(a) 



(b) 



(c) 



FIG. 2: Drell-Yan scattering amplitudes via (a) q + q ^ J* + g and (h) q + q + g^>-J* + g; (c) A 
typical diagram, from the interference of the amplitudes in (a) and (b), that gives a contribution 
to the SSA. 

when the lepton pair is produced in the forward direction of the polarized beam (at large 
rapidity), so that large momentum fractions in the polarized proton are probed. 

As we discussed in the Introduction, we are primarily interested here in the Drell-Yan 
pair produced at large transverse momentum q± ~ Q. For this to happen, a quark or gluon 
jet must be produced in the hard partonic process against which the Drell-Yan pair recoils. 
In the unpolarized case, to lowest order in perturbation theory, there are two partonic 
subprocesses of this kind: the Compton-type quark-gluon scattering g + (? — > 7* + g, and 
quark-antiquark annihilation q + q ^ 'y* + g. One therefore arrives at the following lowest- 
order (LO) expression of the unpolarized Drell-Yan cross section for producing a virtual 
photon of invariant mass Q, rapidity y, and transverse momentum q±: 



d^a 



dQ'^dyd?qx 
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X(5(s + t + u — Q 



(5) 



where a^ is the strong coupling constant, and ctq = 4,'na1^J?>NcsQ'^ with the electromagnetic 
coupling aem and the hadronic center-of-mass energy squared s = {P + P')"^. Furthermore, 
X and x' are the partons' momentum fractions, and the partonic Mandelstam variables are 



defined as s = {xP + x'P'Y, i = {xP — q^* 



u 



[x'P' — q^*y with the virtual photon 



momentum q^* . Note that we have omitted the scale dependence of the parton distributions, 
as we will do throughout this paper. 

In order to have a non- vanishing single-transverse-spin asymmetry in a hadronic process. 



an interference of two amplitudes with different strong-interaction phases is required. When 
the underlying scattering mechanism is hard, as is the case for our Drell-Yan observable at 
large Q and g^, the difference in strong interaction phase may arise from the interference 
between a real part of the scattering amplitude in Fig. 2(a) and an imaginary part of the 
partonic scattering amplitude with an extra gluon in Fig. 2(b). The interference of these two 
amplitudes corresponds to the typical diagram contributing to the SSA shown in Fig. 2(c). 
The additional gluon of momentum kg from the polarized proton, which leads to the twist- 
three quark-gluon correlation function as shown in Fig. 1, can attach to any propagator of 
the hard part represented by a light circle in Fig. 2(b) or (c). The imaginary part of the 



amplitude with an extra gluon is provided by the pole of the parton 



propagator associated 



Q 



In the previously 



with the integration of the gluon momentum fraction Xg = k'^ / P^ 
considered cases of the SSAs in direct-photon or inclusive-hadron production, such poles only 
occur when the additional polarized gluon has a vanishing momentum, x^ = (xi = X2), 
while it is attached to an external on-shell parton 71. These poles were referred to as "soft 

nn 

poles" |2^|2a]. For example, if the polarized gluon attaches to the incoming unpolarized anti- 
quark, the on-shell propagation of the anti-quark line will generate such a soft and un-pinched 
gluonic pole. However, for the Drell-Yan process, there are two observed hard scales, Q and 
g_L, because of the outgoing photon being off-shell, and there could be additional imaginary 
contributions associated with poles of propagators in the partonic scattering amplitude for 
which Xc, 7^ (xi 7^ X2). Such poles were called "hard poles" J2^, |2a|. In the calculations 
presented below, we shall consistently include both contributions, soft-pole and hard-pole 
ones. Compared to the unpolarized case in Eq. ©, we now find the following structure for 
the LO single-spin cross section: 

wi ^ -"'"^-'-li Jt'^T.< m + Hi) «(-') + {HI + H^,) m 

x5{s + t + u-Q'^) , (6) 

where H^l and H^ are the soft- and hard-pole contributions for scattering off an unpolarized 
anti-quark or gluon, respectively. Each of these will be a function of the partonic Mandelstam 
variables, of x and x', and of the twist-three quark-gluon correlation function Tp(xi,X2). 
Note that x is the same as in the unpolarized case, determined in terms of the external 
hadronic variables through the delta- function in ®. As will be described in detail below, 
Xi and X2 depend on x and are fixed by the on-shell conditions set variously by the soft and 
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FIG. 3: A typical Feynman diagram (a) and its complex conjugate (h) contributing to the single-spin 
asymmetry for the Drell-Yan process through the quark- antiquark scattering channel, containing 
both soft-pole "(S)" and hard-pole "(H)" contributions. The bars indicate the places where a pole 
arises and generates a strong interaction phase. 

hard poles, and by that for the unobserved final-state parton in the hard process. 

In the following subsections, we will first discuss the generic features of the soft and hard 
poles, and then present their respective contributions to the single-transverse polarized cross 
sections for the quark-antiquark and quark-gluon scattering channels. 

A. Generic Structure of Soft and Hard Poles 



We will take the diagrams shown in Fig. 3 as examples of how the soft and hard poles 
arise. Let p' denote the momentum of the incident anti-quark, kqi that of the initial quark to 
the left of the cut, kg2 that on the right, and kg = kg2 — kgi the momentum of the polarized 
gluon attaching to the hard part. This attachment may happen on the left side of the cut, 
as shown in Fig 3(a), or on the right side as in 3(b). In order to analyze the pole structure 
in the scattering amplitudes, it is convenient to consider only the dominant components of 
the relevant momenta. For example, p' may be chosen to only have a light-cone "minus" 
component, while kgi and kg2 are then dominated by their plus components. When the 
polarized gluon attaches to the left side of the cut in Fig. 3(a), we find that the on-shell 
condition for the gluon radiated into the final state fixes k'^2 (o^^ equivalently, X2) to be 
X2 ~ X, while kt'i {xi) is determined by the (soft- or hard-) pole condition. On the other 
hand, the momentum fiow in Fig. 3(b) is opposite: Xi is determined by the on-shell condition 
for the outgoing gluon, and X2 by the pole of a propagator. 

In Fig. 3, we have denoted the poles by bars, and introduced the labels (5) and (H) for 



the soft and hard poles, respectively. A soft pole arises in Fig. 3(a) on the anti-quark line 
carrying momentum p' + kg, because we have 

I = -J: L. (7) 

{p' + kgY + ie 2p'-P+ Xg + ie ' ^ ^ 

which provides a phase proportional to S{xg) = 6{xi — X2)- Therefore, the quark-gluon 
correlation function Tf{xi,X2) will be probed when both of its arguments are identical. 
After combining with the on-shell condition for the unobserved final-state gluon, one finds 
in fact that the contribution enters with Tf{x, x) ^. Similarly, the soft pole in Fig. 3(b) will 
arise through 

' ' (8) 



(p' — kgY — i^ 2p' P+ —Xg — ie 
so that it enters with opposite sign compared to that of Fig. 3(a), but with otherwise the 
same quark-gluon correlation function. Therefore, the contributions by the two poles in 
Eqs. dZI),® may be combined. 

Next, we turn to the hard pole in Fig. 3(a). It occurs when the quark propagator carrying 
momentum kqi — q^* goes on mass-shell. The propagator reads 

I ^ I (9) 

{kgi - q^*f + It -2A;+g- + Q^ + le ' ^ ' 

where we have neglected the transverse momentum fcgi^. It is easy to see that this propagator 
has a pole at k'^^ = xP^Q^/{Q^ — t), corresponding to Xg = Xg, where 

Xg = -xi/{Q^ - i) . (10) 

Because t = {xP — q-y^Y < 0, this pole is not forbidden kinematically. Furthermore, Xg > 0, 
which is the reason why we refer to the pole as "hard". After taking the pole, we find 
Xi = X — Xg. Using the on-shell condition for the outgoing gluon, the associated twist-three 
correlation eventually is Tf{x — Xg,x), different from that for the soft-pole contributions. 
Likewise, the hard pole in Fig. 3(b) arises from the propagator 

{kg2-qr)^~te "^ -2k^2<ir + Q^ ' ^^ ' ^^^^ 

and is at k^2 = xP^Q^/{Q^ —t), corresponding to Xg = —Xg = xt/{Q'^ — t) < 0. The 
associated twist-three quark-gluon correlation function becomes Tp{x,x — Xg). Because of 



^ As we shall discuss later, there are also contributions with the derivative dTp{x,x)/dx. 
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the symmetry property (JH) of the function, we have Tf{x, x — Xg) = Tp{x — Xg, x), and we 
may again combine the contributions from the two sides of the cut in the diagram. 

From kinematics, we find that hard poles only arise from t-channel propagators, and only 
on the side of the diagram that contains the additional gluon from the polarized proton. 
For each diagram, we have to include all contributions by the soft and hard poles, in order 
to obtain the complete result. Although we will calculate their contributions separately, the 
soft and hard poles do overlap in some kinematical regions, for example, when the transverse 
momentum q± is much smaller than Q. We will discuss this more in Sec. IV. 

B. Soft-Pole Contributions 

In this subsection, we calculate the soft gluonic pole contributions to the single- 
transversely-polarized Drell-Yan cross section in the LO quark-antiquark and quark-gluon 
scattering channels. The calculation is rather similar to that performed for direct-photon 
production in Ref. [3], except for the fact that the photon is now virtual. 

1. Soft Poles in Quark- Antiquark Scattering 

We start by recalling the unpolarized hard-scattering function for the subprocess q + q ^ 
•J* + g, appearing in Eq. (0): 

aUsXu) = 2cJi + - + %^) , (12) 

\t u ut J 

where Cp is the color-factor and the partonic Mandelstam variables are as defined after 
Eq. (0) at the beginning of this section. 

For the single-spin cross section in this channel, there are eight partonic diagrams pos- 
sessing soft poles. Four of them are shown in Fig. 4, and the other four are obtained by 
attaching the polarized gluon in the same way on the right side of the cut. Other attach- 
ments of the polarized gluon either do not give a soft pole, or simply cancel after summing 
over cuts. 

It is useful to recapitulate some of the key steps in calculating the diagrams in Fig. 4 [7| . 
We perform our calculations using a covariant gauge. The polarized gluon is associated with 
a gauge potential A^, and one of the leading contributions comes from its component A^ . 
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FIG. 4: Feynman diagrams for the soft-pole contributions to the single-transverse-spin asymmetry 
through the quark- antiquark scattering channel. The bars indicate the propagators yielding poles. 



The gluon's momentum is dominated by XgP -\-kg±, where Xg is the longitudinal momentum 
fraction with respect to the polarized proton. The transverse momentum kg± flows through 
the perturbative diagram and returns to the polarized proton through the quark lines. The 
contribution to the single-transverse-spin asymmetry arises from terms linear in kg± which, 
when combined with A^, yield d-^A'^, a part of the gauge field strength tensor F-^~^. In 
order to compute this contribution, we expand the partonic scattering amplitudes in terms 
of kg± up to the linear term. One important contribution to the kg± expansion comes from 
the on-shell condition for the outgoing unobserved quark, whose momentum depends on 
kg±. As was shown in |i2|, this leads to a term involving the derivative of the correlation 
function Tp. Apart from this contribution, the kg± expansion of the propagators and the 
Dirac traces yields a term proportional to the correlation function itself. This is the so-called 
non-derivative term. 

The sum of the soft-pole contributions by the diagrams in Fig. 4 and their "mirror images" 
gives the function H^ in Eq. (jHl). We find: 



H', 



x-—Tf{x,x) 



_J!1+Tf{x,x)^^ 
—u —u 



(13) 
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FIG. 5: Feynman diagrams for soft-pole contributions to the single-transverse-spin asymmetry 
through quark-gluon scattering. 



where the hard coefficients D^^ and A^^^ are 



D 



m 






2{Nl - I 



■Crqq{s,t,u) 



2Nc t^u 



[Q\u^-P) + 2Q^s{Q^-2i) 



u 



m , 



where dqq in the ffist hne is the unpolarized hard-scattering function given in Eq. ()12j) . 
implying that the hard-scattering function for the derivative term differs from the unpo- 
larized one only by the color factor 27[. In the real-photon limit, Q^ — > 0, we obtain the 
annihilation contribution to the direct-photon single-spin cross section. 



2. Soft Poles in Quark- Gluon Scattering 

Again we start by giving the hard-scattering function for the unpolarized case: 

2Tr 



aqg{s,t,u) 



s -t 2Q^u 
-i s si 



(14) 



where Tr = 1/2. 

Feynman diagrams for the soft-pole contributions to the single-transverse-spin asymmetry 
are shown in Fig. 5, where again the diagrams for the attachments of the polarized gluon 
to the incident gluon on the right side of cut have been omitted. Following the same 
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calculational procedure as outlined earlier, the soft-pole contributions for the quark-gluon 
channel are: 



H', 



x—Tf{x,x) 



'^^ +Tp{x,x)^ , (15) 



—u —u 

where the first term is the derivative term and the second the non-derivative one, with 

iv; = 2^Ni-i) ¥-s f^'^^' ' ''^ + ^^'^^^' " ^^ ~ ^^ + ''^'] • 

We note that apart from the color factor, the hard coefficients D'^ and A*"* for the quark- 
gluon scattering channel can be obtained from those for the quark-antiquark channel by 
"crossing" s ^^ u. In the limit Q^ = 0, we obtain the direct-photon cross section, which 
was calculated previously in [7J. We agree with the derivative term found there, but find a 
difference for the non-derivative one. 

C. Hard-Pole Contributions 

In this subsection, we will calculate the hard-pole contributions to the single-spin cross 
section. The calculational procedure is similar to that for the soft-pole contributions in the 
previous subsection. We will again discuss the quark-antiquark and quark-gluon channels 
separately. 

1. Hard Poles in Quark- Antiquark Scattering 

In Fig. 6 we show the hard-pole diagrams for the quark-antiquark scattering channel, 
again omitting the "mirror" diagrams. As discussed in subsection III. A, only t-channel 
propagators give a hard pole. In the quark-antiquark scattering channel, there are only 
two basic diagrams containing a t-channel propagator, shown in Figs. 6(a) and (d). How- 
ever, because of gauge invariance, we need to allow all possible gluon attachments in these 
diagrams. We thus have attachments to the initial antiquark line, to the outgoing gluon 
line and to the t-channel quark propagator, as shown in the figure. Therefore, including 
the diagrams with gluon attachments on the right side of the cut, there will be a total of 
twelve diagrams containing hard poles in the quark-antiquark scattering channel. We may 
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FIG. 6: Hard-pole contributions in the quark- antiquark scattering channel. The bar indicates the 
propagator that produces the pole. 

group their contributions into two color structures: Tr (T"T''T'^T*) for diagrams (a) and 
(d), and ^^{T^'T'^T^T^) for diagrams (c) and (f). The first color factor is the same as for 
the soft gluonic pole diagrams. The color traces for diagrams (b) and (e) can be decom- 
posed into these two color structures, and combined with (a,d) and (c,f), using the fact that 
ifabcT^ = T^T^ - T^T". 

The expansion in kg± for the hard-pole contributions is somewhat more complicated 
than for the soft-pole ones, because the position of the hard pole itself depends on kg±. 
As an example we consider the diagram in Fig. 6(a), for which a hard pole arises from the 
propagator 

{kqi - q-y*y + ie -2k'^^q~, + Q^ - 2kqij^ ■ q^,± + ie ' 
where we have kept the full dependence on the transverse momentum. In the fc^x-expansion, 
the kqi± term in the above denominator will make a contribution to the term linear in kg±, 
because kg± = kg2± — kgi±. This will lead to a derivative term (of Tp), similar to the double- 
pole contributions found in [7] for the SSA in single-inclusive hadron production. In the 
calculation of this derivative term, we may set kg± = everywhere else in the scattering 
amplitude. Because all the propagators with hard poles in Fig. 6 are t-channels and hence 
have the same momentum, namely kgi — q^* , the expansion will be the same for all diagrams. 
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and we may add their contributions. After summing over all the diagrams, we find that the 
total derivative contribution vanishes. More precisely, the derivative contributions from 
the hard poles in diagrams (a), (b), (d) and (e) (which all contain the color structure 
Tr(T"T^T"T^)) cancel each other, and the same happens for the contributions by diagrams 
(b), (c), (e) and (f) (which enter with Tr(T"T'*T^T^)). 

The expansion of the delta function associated with the on-shell condition for the outgoing 
unobserved gluon {{kq2 + p' — q-y*)^ = 0) also contributes to the derivative term for each 
individual diagram. However, after summing over all diagrams, again the net derivative 
contribution vanishes. In summary, we do not have any derivative contributions originating 
from the hard poles. 

For the non-derivative contributions, we need to do an expansion of the numerator (for 
example, the Dirac trace) and of the other propagators in the squared amplitude. Apart from 
that, the derivative terms just discussed will make contributions to the non-derivative pieces. 
Our way of dealing with the non-derivative contributions is to first evaluate the scattering 
amplitudes with their full dependence on /c^^, and then to use the on-shell condition for 
the outgoing unobserved gluon and the hard-pole condition, in order to fix the light-cone 
momentum fractions of the quarks including their dependence on /c^^. For example, the 
on-shell condition for kq2 + p' — k^* in Fig. 6(a) determines the value for kt2: 

'2kq2±.q±^ 
u 

Similarly, the hard-pole condition {kqi — A;^.)^ = fixes the value for kq^ as 

From the above equations, we see that both k'^^ and /ci depend on /c^_|_, so that their 
expansion will contribute to a non-derivative term. 

Summing over all the contributions by the diagrams in Fig. 6 and their "mirror im- 
ages" , we obtain the final result for the non-derivative hard-pole contribution for the quark- 
antiquark channel: 

^2t.2 r T si 

(19) 



k-=xP-{l---^^^^^^\ . (17) 



Hq = Tf{x - Xg, x) X ^ ^^ 



' +c,- ' 



2Nc s + u 



where as before Xg = —xi/{Q'^ — i) is different from zero, reflecting the hard-pole nature of 
the contribution. The two color factors in the above equation are associated with the two 
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FIG. 7: Hard pole contributions in the quark-gluon scattering channel. The bar indicates the 
propagator that produces the pole. 

independent color structures discussed earlier, the first one coming from Xr(T'*T^T"T^) and 
the second one from Tr (T^T^T^T''). 

2. Hard Poles in the Quark- Gluon Channel 

By a similar procedure, we calculate the hard-pole contributions for the quark-gluon 
scattering channel. In Fig. 7, we show the relevant diagrams. Following the same arguments 
as given for the quark-antiquark scattering channel, only t-channel diagrams contribute to 
the hard poles, and their total contribution is 



H^ = Tf{x - Xg,x) X 



—ut'^s 



-Nh 



2(iV< 



+ T, 



' s + u 



(20) 



Again, there are no derivative terms, and the hard coefficient for the non-derivative term is 
related to that for quark-antiquark scattering given above by s ^^ m crossing. As for the qq 
case, there are two contributing color structures, one associated with ifabcT^T^{T"'T''T'^) (for 
diagrams Fig. 7(a), (b), (d), and (e)), and one with 'Xr(T"T'^T^T^) (for diagrams Fig. 7(c), 
(b), (f), and (e)). 
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IV. LOW TRANSVERSE MOMENTUM LIMIT 



Having calculated the full LO single-spin cross section for the Drell-Yan process at large q± 
within the ETQS approach, we are now in the position to investigate the limit Aqcd ^ Q'± *^ 
Q, trying to make contact with the approach based on transverse- momentum dependent 
factorization. To this end, we expand the formulas given in the preceding section for small 
Ql/Qi extracting their leading contributions. 

The partonic Mandelstam variables may be written as 



<il 



:i-6)(i-6) 



u 



1-6 



(21) 
(22) 
(23) 



1-ei ' 

where ^ = zi/x and ^ = Z2/x' , with zi = Q/y/se^ and Z2 = Qj \fse~y . The phase-space 
delta-function for the on-shell condition for the outgoing parton takes the form 



where the "plus" -prescription is defined in the standard way through 



(24) 



\z- ^(^) 



'rf^M_/(l)+/(i)in(l-x), 



(25) 



;i-^)+ 

for any suitably regular function /. 

After applying Eq. (jUj) to the results presented in the previous subsections, we find for 
the small-g^ behavior of the unpolarized cross section for q + q —>■ '~f*g: 



dQ'^dyd'^q± 



a 



1 C dx dx' 



/vr"^ q^ J XX' 



rX]^?^'^^)^^^') 



;i-e 



-m - 1) 



i;+ 



(1-6;+ 

For q + g ^ ■y*q scattering, we obtain 



91 J 



(26) 



dQ'^dyd'^q± 



a., 1 



CTO 



J^i^-^J24<l(^)9ix')6il~^,)[a + {l-^2n . (27) 



47r2 qj_ J XX' 



These results are well-known in the literature 



Q. 
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We now make similar expansions for the single-spin cross sections. As we anticipated, 
at low transverse momentum, the soft- and hard-pole contributions overlap: at q± = 0, the 
twist-three quark-gluon correlation function associated with the hard-pole contributions, 
Tp{x — Xg,x), becomes Tp{x — x{l — ^i) / ^2, x) , which at ^1 = 1 is identical to the correlation 
function Tp{x,x) that accompanies the soft poles. 

For the qq channel, we find: 



^4^^99^7*9 (5j 



dQ'^dydPqi 



Oq e°''^S±a 



#^/f^*''(*(«-^'-^-^*«'-^w-(^«) 



where 



A 



1 



2N, 



c 



x§^Tp{x,x) 



[1 + ^l) +Tf{x,x -Xg 



1 + ^1 



ll)+ 



B 



-|-Tp(x, x)- 
CfTf(x,x 



1-6)^(2^1 + 1) -2 
(1-6)+ 



\ + CfTf{x, X - Xg)- -\ 



:i - 6)- 



nin 






(29) 
(30) 



with Xg = {1 — C,i)x. 

For the qg channel, the contribution involving the derivative dTF{x,x)/dx turns out to 
be of higher order in q±/Q. The full non-derivative term becomes 

d^Aa^a^^'^iS^) 



dQ'^dyd'^q± 



aoe"^^Xa7%$^e;T^(^i,^i) 



W. 



x-^tW— ^(x2)fe^ + (l-6)^] 
X2 



2tt^' 



(31) 



In the following section, we will reproduce the above results from a QCD factorization in 
terms of transverse-momentum-dependent parton distributions. 

V. TRANSVERSE-MOMENTUM-DEPENDENT QCD FACTORIZATION 



When q± ^ Q, the Drell-Yan cross section in leading order in q±/Q may be calculated 
from a QCD factorization theorem involving transverse-momentum-dependent (TMD) par- 



ton distributions 



m 



Q. 



Therefore, one expects that the results given in the previous 



section can be reproduced by this alternative approach. Indeed, for the unpolarized cross 
section this is well established Il3l ^ . 
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A. Transverse-Momentum-Dependent Parton Distributions 



Transverse-momentum-dependent parton distributions were first introduced by Collins 
and Soper [ij]. They provide more information about the structure of the nucleon than 
is contained in the usual Feynman parton distributions. Based on rigorous factorization 
theorems [ij, llJ, ll5|) these distributions may be extracted from various processes that 
are characterized by a large momentum scale and a much smaller, measured, transverse 



momentum, such as semi-inclusive DIS and the Drell-Yan process. In this paper, we 



bllow 



a definition of the TMD distributions in Feynman gauge with explicit gauge links [17[. We 
avoid the axial gauge because of the complications presented by additional gauge links at 
space-time infinity [l8|. 

The TMD quark distributions of a polarized proton may be defined through the following 
matrix: 



M"^ = P+ 



dC 



-ix(-p^ 



Jj^^b^-k^ 



PS 



<(r,o,6x)>^."(o) 



PS 



(32) 



277 ' J (27r)2 

where the vector P = {P^, 0~, 0^) is along the momentum direction of the proton, S is the 



polarization vector, and ^v{0 i^ defined as 



^,(0 = /:.(-oo;0V^(0 



(33) 



with the gauge link 



£„(-oo;0 = exp 



-w 



d\v ■ A{\v + ^) 



(34) 



This gauge link goes to — oo, indicating that we adopt the definition for the TMD quark 



distributions for the Drell-Yan process |l6lll7lll8||. In a coordinate system where P+ ^ P^, 
the vector v in the above equations is taken to have v~ ^ v~^, v± = 0, with a nonzero 
f + component that regulates the light-cone singularities. The leading-power expansion of 
the density matrix 7W contains eight quark distributions [10|. Here, we are only interested 
in the TMD distributions for (un)polarized quarks and antiquarks in the polarized proton, 
in particular in the Sivers functions. The physics of the Sivers functions may be viewed as 
follows. Consider a transversely polarized proton with large momentum P. The distribution 
of quarks with longitudinal and transverse momenta xP and k± can have a dependence on 
the orientation of k± relative to the proton's polarization vector S±. Keeping only the 
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FIG. 8: Feynman diagraTns contributing to the spin-independent quark distribution at large trans- 
verse momentum k± ^ AqcDj arising from the integrated quark (a-d) and gluon (e) distributions. 
Double lines denote Eikonal lines. 

unpolarized quark distribution and the Sivers function, we have the following expansion for 
the matrix M.: 



M^- 



IVIp 



(35) 



where q{x, k±) is the TMD distribution in an unpolarized proton, qrix, k^) is the Sivers 
function, and Mp is a hadron mass, used to normalize q{x,k±) and qT{x,k±) to the same 
mass dimension. 

When the transverse momentum k± is low, say, of order of Aqcd? the TMD parton dis- 
tributions are entirely non-perturbative objects. Thus, the only way of treating them is to 
parameterize them in a suitable way, and to fit them to data. However, for much larger 
transverse momentum, k± ^ Aqcdj the fc^-dependence of the TMD parton distributions 
can be calculated within perturbative QCD, because a hard gluon has to be radiated in order 
to generate the k±. Through the pQCD calculation, the TMD parton distributions may be 
related to the fc^-integrated parton distributions (or, in case of the Sivers functions, to the 
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FIG. 9: Feynman diagrams contributing to the Sivers functions at a large transverse momentum 
k± ^ AqcD; arising from the twist-three quark-gluon correlation in Fig. 1: soft-pole contributions. 

twist-three quark-gluon correlation function, as we shall show below), multiplied by calcu- 
lable hard coefficients. For example, the TMD quark distribution for an unpolarized proton 
at large fc^ will receive contributions from the fc^-integrated quark and gluon distributions. 
The relevant Feynman diagrams are shown in Fig. 8, and the results for the first-order term 
are (see, for example, |l4l3fll|): 

1 ^ C dx , . 



q{z,k±) 



a« 



2n^kl 
as 

27T^'k 



C 



-q[x) 



1-0 



^^^%.(e-l)fln.^^^^ 



k' 



(36) 



where q{x) and g{x) are the integrated quark and gluon distributions, ^ = z/x, and ("^ = 
{2v ■ PY /v"^. A similar expression is obtained for the anti-quark TMD distribution. We note 
that upon including virtual corrections, we would find that the term l/k^ would be converted 
to a term l/{kY)+, and similarly for the logarithmic term. This "plus" -prescription would 
make the distribution integrable at low fc^; for details, see for example |8l| . 

In the same spirit, the Sivers function at large k±^ can also be calculated perturbatively. 
Because it is (naively) time-reversal-odd, the only contribution comes from the twist-three 
quark-gluon correlation, with the phase provided by the hard part. This is very similar to 
the case of the single transverse-spin dependent Drell-Yan cross section calculated in Sec. III. 
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Also in the present calculation, we will have soft-pole and hard-pole contributions, shown 
in Figs. 9 and 10, respectively. Carrying out the calculations accordingly, we find for the 
soft-pole contributions: 



qriz.k^ 



I soft— pole 



as 2Mp 1 



+ Tf(x,x) 



dx 

X 

2MC , 1^ 2 



^-q^Tf{x,x) 



(1 + a 



:i-o'(2^+i) 



;i-o 



+ 



5(e-l)(ln^-l 



(37) 



where the 1/(A;^)^ behavior also follows from power counting, and where ^ and ( as defined 
above after Eq. (l36l) . We note that our choice of the direction of the gauge link in Eq. (l3Hl 

nrfn 

is crucial since it determines the sign of the Sivers function |lfil . I17L I18| . Only with the 
correct choice (in this case, to — cxd) can the factorization work out eventually. In case of 
semi-inclusive DIS a different choice is necessary, and the resulting Sivers function will have 
opposite sign. 

The hard-pole contributions can be calculated similarly, and we find 



qTiz,k±)\ 



hard— pole 



as 2Mp 
1 + e 



Cf + 



2Nc, 



X 



ln3!-l 



(38) 



where Xg = (1 — ^)x. Adding the above soft-pole and hard-pole contributions, we obtain the 
final result for the Sivers function at large k^: 



'-(^'•*J = i^f;^/T{-^ + ^-^-<^-^'*'«'-')('"ir-' 



(39) 



where A has been given in Eq. (jSH))- 

B. Unpolarized Cross Section 

The spin-independent cross section for the Drell-Yan process has the following TMD 
factorization in the regime q_L <^ Q [iJ, llJ| : 
d'^a 



dQ'^dyd'^q± 



(To ^ ej / d%^d^k2±d'^X±6^^\ki± + hx + X± - q±) 



X g(zi, A;i^, Ci)g(^2, k2±, (2) H{Q'') {S{\^, p)) 



-1 



(40) 



where as before q{z, k±) and q{z, kj_) are the TMD quark and antiquark distributions. H is 
a hard factor and is entirely perturbative. It is a function of Q ^ g_L only. The soft-factor S 
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FIG. 10: Same as Fig. 9, hut for the hard-pole contributions. 

is a vacuum matrix element of Wilson lines and captures the effects of soft gluon radiation. 
Since the soft-gluon contributions in the TMD distributions have not been subtracted, the 
soft factor enters with inverse power. As mentioned earlier, in order to regulate the light-cone 
singularities, we introduce the off-light-cone vectors Vi and V2 for the TMD distributions. 
We also define Ci = (2f 1 ■ P^/vf and Q = {'^'^2 ■ P'Y /vl, and the soft-gluon rapidity cut-off 



p = A/(2t'i ■ V2Y /v'lv2. In a special coordinate frame, zfCl = Z2C2 = PQ'^ [14]. There are 
also explicit renormalization scale dependences of the various factors in the factorization 
formula, which have been omitted for simplicity. 

In Eq. flHfij) in the previous subsection, we have already given the unpolarized TMD 
quark distribution, calculated perturbatively for kj_ ^ Aqcr^Also the soft factor may be 
calculated in this fashion. To leading order in a^, one finds 14j |. 

1 



S{X±) = 6^'\\^) + 



tt.9 



27r2A^ 



3, 

-Cf (lnp2 - 2) 



(41) 



By adding virtual contributions, the soft function is normalized in such a way that it obeys 
Jd^X±S{X^) = 1. 

Inserting the perturbative TMD distribution (pH|l and the soft function ()4H1 into the fac- 
torization formula (001), "we find that indeed the unpolarized cross section given by Eqs. ()27p. 
fl26|) is reproduced, including the quark-gluon scattering piece. Here we use the above nor- 
malization condition for the soft function, and the normalization that the integration over 
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the TMD distribution yields the normal Feynman parton distributions, 

(fk±q{zi, k±) = q{zi) , / (fk±q{z2, k±) = q{z2) ■ (42) 



C. Polarized Cross Section 

In |lj] also a factorization formula for the transverse-spin dependent Drell-Yan cross 
section aX q± <ti Q was established, which reads: 

dQ^dyd^q± Mp J q]_ 

xqT{z,,k,^, Ci) qiz2,k2^X2) HiQ') (SiX^))-' , (43) 

where qT{z,k±) is the Sivers function defined in Eq. (jH3|) . q{z,k^) is again the anti-quark 
TMD for the unpolarized proton, and H and S are the hard and soft factors introduced 
earlier. 

We have eriven the first-order perturbative result for the Sivers function in Eq. (1391) . As 

n nn 

was shown in |19| (see also |22,|23|), its A;5.-i^oi^6nt is related to the twist-three quark-gluon 
correlation function defined in Eq. Q of Sec. II: 

— / d'^k±_k'iqT{x,k±) = Tf{x,x) . (44) 

For deep-inelastic scattering, this relation remains true, except that the sign will be opposite. 
Inserting finally all TMD functions into the factorization formula ()43|1 . we reproduce the 
differential Drell-Yan single-transverse-spin cross section given at low transverse momentum 
in Eqs. (jS^, ()28p of Sec. IV. This concludes our demonstration that the ETQS mechanism 
and the approach based on TMD factorization produce identical results in the kinematic 
regime where they overlap. 

VI. SUMMARY AND OUTLOOK 

In summary, we have studied the single-transverse-spin asymmetry in Drell-Yan pair 
production at large and moderate transverse momenta. At large transverse momenta qj_ ~ 
Q, the asymmetry is power-suppressed by 1/Q. Therefore, the coUinear-factorized approach 
based on twist-three quark-gluon correlations should provide the appropriate description. 
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We have used this approach to derive the single-spin asymmetry. We have then expanded 
the result for moderate transverse momenta of the lepton pair, Aqcd ^ q^ ^ Q- In this 
regime, one knows that a factorization theorem in terms of transverse-momentum dependent 
parton distributions holds, involving in particular the Sivers function. We have verified that 
indeed a smooth transition from the higher-twist mechanism to the one based on TMD 
factorization occurs, with the two approaches describing the same physics in the region 
of overlap. This unifies the two approaches, and provides input for new strategies in the 
analysis of experimental data for single-spin asymmetries. 

In p], it was noted that there are additional twist-three quark-gluon correlation functions 
that may contribute to the SSAs in hard processes. One of them is a spin-averaged and 
chiral-odd correlation, defined by [7[ 

tP{x,,X2) = y ^^L^e^-i^+^-e^(-2-x,)P+C- (p|^(o)a+"^F+"(C-)^(r)|P> . (45) 

When combined with the quark transversity distribution for the polarized proton, the above 
correlation in the unpolarized proton can also contribute to the SSA for the Drell-Yan process 
at large transverse momentum. However, after a straightforward calculation we find that 
the contribution vanishes in the low transverse momentum limit q± <C Q. This observation 
is consistent with the TMD factorization at low transverse momentum [lj|, where only the 
Sivers function contributes to the SSA for the Drell-Yan process after integration over the 
lepton angles. 

There are a number of directions in which our work could be extended. First, an in- 
teresting topic will be to study the single-spin asymmetry in semi-inclusive deep inelastic 
scattering (SIDIS) ep -^ eirX [32], again at large and moderate transverse momenta of the 
final-state hadron. Also here, a connection between the twist-three mechanism and the TMD 
factorization formalism can be established. The Sivers function calculated perturbatively 
in this paper may be immediately applied to SIDIS, after changing its sign. Although the 
calculation for SIDIS should be straightforward following our derivations for the Drell-Yan 
process in this paper, the details will differ. It will be crucial to demonstrate that the two 
approaches are indeed also consistent for the SIDIS process at Aqcd "^ q± "^ Q, and to 
explicitly verify the process-dependence for the single-transverse spin asymmetries in hard 
processes. It also needs to be demonstrated that the Collins function mechanism J22| can be 
reproduced in the low-g^ region. We will present a study of SIDIS in a future publication. 
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Second, the consistency between the twist-three and the TMD factorization approaches 
should also lead the way to resumming the large logarithms from higher-order soft gluon 
radiation. As shown in Eq. (j28|) - (j3(J|) . the first-order expressions contain a large logarithm 
oc \iiQ'^/q\. This large logarithm will be present in higher orders, with two additional 
powers of the logarithm at every new order in perturbation. In order to have a reliable 
theoretical calculation, these logarithms need to be resummed. The resummation procedure 
may follow the classic example in [l3| for the unpolarized Drell-Yan cross section. The 
relevant CoUins-Soper equation for the energy evolution of the spin- dependent TMD quark 
distributions has already been derived in J2J|. The resummation effects will be particularly 
important at small q^jQ- 
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